Visibility of ultracold Bose system in triangular optical lattices 
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In this paper, by treating the hopping parameter in Bose-Hubbard model as a perturbation, with 
the help of the re-summed Green's function method and cumulants expansion, the momentum distri- 
bution function of the ultra-cold Bose system in triangular optical lattice is calculated analytically. 
By utilizing it, the time-of-flight absorption picture is plotted and the corresponding visibility is 
determined. The comparison between our analytical results and the experimental data from Ref.pJ 
exhibits a qualitative agreement. 
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I. INTRODUCTION 

During the past decade, the properties of ultra-cold atomic gases in optical lattices jilj have received a great deal of 
attention on account of their novelty and various potential applications [2|]. In the case of bosons, the delicate balance 
between the atom-atom on-site repulsion and the hopping parameter leads to a Mott insulator-Superfluid (MI-SF) 
quantum phase transition which has been confirmed experimentally Due to the precisely controllable parameters, 
the optical-lattice-ultra-cold-atom systems play roles of test ground for quantum many-body systems in condensed 
matter physics, this is the so-called quantum simulation Q. 

Beside the simple cubic optical lattices, in which almost all experiments with ultra-cold atoms to date have been 
performed due to the ease of the ease of the implementation, the triangular and hexagonal optical lattices have also 
been recognized recently In fact, due to the complexity of the lattice structure in these systems, novel and rich 
new phases will be exhibited, for instance, the effects of geometrical frustration have been achieved in triangular 
optical lattice cold atomic system experimentally [!, H|. Hence, it is worth to investigate these systems, especially 
their quantum phase transitions, in a systematic way. 

In our previous work [7], we have calculated the phase boundaries of MI-SF quantum phase transitions of ultra-cold 
bosons in triangular, hexagonal, as well as Kagome optical lattices analytically via the method of the field-theoretical 
effective potential [8[ , the relative deviation of our analytical results from the numerical results [|] is less than 10%. 
However, these results cannot be compared with the experimental observation, since the phase boundaries of the 
ultra-cold system are not able to be detected directly in experiments. Instead, time-of-flight measurement @], which 
reveals the momentum distribution of the system, is the standard experimental technique for investigating cold atom 
systems in optical lattices [HI [H| • 

In this paper, by treating the hopping parameter as a perturbation, with the help of Green's function calculation (l2l . 
Il3j , we are going to calculate the time-of-flight absorbtion pictures and the associated visibility as well analytically for 
ultra-cold scalar bosons in triangular optical lattice. Our results exhibit a qualitative agreement with the experimental 
data. 

II. THE MODEL 

As is well known, a system of scalar bosons trapped in a homogeneous optical lattice is portrayed by the many-body 
Hamiltonian which takes the form as 

H = J dr^(r) (-t^V 2 + V(r)j if,(r) + | J dr^ (r)^ t (r)^(r)V>(r), (1) 

with ip(r) (^>t(r)) being annihilation (creation) bosonic field operators and V(r) being the optical lattice trapping 
potential. The parameter g is the interaction strength between two atomic particles. 
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In fact, due to the modulation of the optical lattice potential wells, in the extremely low temperature limit, a 
single band approximation is adequate, the bosonic field operators i/^( r ) an d "0( r ) can be expanded in the basis of 
orthonormal Wannier functions w(r — r 3 ) of the lowest band [H], HH as 



3 

^ + 0) = ^w*{T-Yi)a\, 



(2) 



where a\ (di) is the creation (annihilation) operator of scalar boson on site of i. 

By utilizing the above expressions, the Hamiltonian in Eq.([T]) is reduced in the tight-binding limit to the Bose- 
Hubbard Hamiltonian 



Hbh — H i + H 
with the diagonal part Hq and the hopping part Hi being 

U 



H o = Y ^iifH - 1 ) - Hi 



(3) 



(4) 



respectively, where n, = a\di is the particle number operator on site i. J is the hopping amplitude for the bosons 
between the nearest neighbor sites i and j while U denotes the on-site repulsion and fi is the chemical potential. As 
we can see from Eq.((T]) to Eq. (J4j , J and U can be calculated explicitly by [lH 



J = — drw*(r — Ti) I — 



2yj2 



2m 



+ V(v) I w(r — rj) 



and 



U = g dr\w(r)\\ 



(5) 



(6) 



respectively, with and Tj being nearest neighboring sites. 

In experiment, the triangular optical lattice is created by three laser beams intersecting in the x-y plane with angle 
of 27r/3 between each other Q, the lattice constant of the obtained triangular lattice is a = 2 A/3, where A is the 
wavelength of the lattice laser beams. The corresponding lattice trapping potential reads 
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Associated to this optical lattice trapping potential, in harmonic approximation, the corresponding Wannier function 
takes form of 



W(T - Ti) 



32V 



exp 



my ^ 

9 } a 2 



((x - Xl f + (y - yi ) 2 ) 



(8) 



with Vq being dimensionless optical lattice depth in unit of recoil energy Er = (h k )/2m, where k = 2tt/X 



III. THE MOMENTUM DISTRIBUTION FUNCTION 

Actually, the quantity that is measured in experiments is the distribution function in the momentum space [l6j 

n(k) = f drdr'e lk '( r - r ')(i/; t (r)V(r')). (9) 
By making use of Eq. ^ , with the help of the Fourier transformation of the operators 

^XXe-* 1 '* a^-^y^e**-'', (10) 
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where N$ is the total of the lattice number, the momentum distribution function is reduced to 

n(k) = N s | iy(k) | 2 (4«k>. (11) 
As is known, the one particle Green's function in the momentum space reads 

G(r|0,k)= <f r [at( r ) k a(o) k ]), (12) 

where r is imaginary time and T T is the imaginary time ordering operator. It is the Fourier transformation of the 
one particle Green's function G(r',j' | r, j) = (T r [a-,(T')a.,(T)]). Apparently, (a k &k) = lim r ioGi(r | 0,k), we 
then translate the problem of calculating the density distribution to a problem of calculating corresponding Green's 
function. 

However, due to the non-commutativity of the two parts of the full Hamiltonian ([3]), the exact information of the 
corresponding eigenstates can hardly be obtained. Nevertheless, the Green's functions can be calculated perturbatively 
via the Dirac representation. To tackle the issue of time-of-flight absorbtion picture which is closely related to the 
MI-SF phase transition, we treat the hopping term in the Hamiltonian (|3|) as the perturbation part, for the opposite 
limit will not lead to a MI-SF phase transition at all (l7j . 

In Dirac picture, the time evolution of operators is only determined by the unperturbed part of the Hamiltonian 
(set Tl = 1) [TH 0(t) = e TH °de~ TH ° , and the corresponding evolution operator takes the form of 

u{/3,0) = T T Uo^^o.i) -rajMajMl . (13) 
Without any difficulty, it can be proved that the one particle Green's function reads 

G(r',f\r,j) = -i 4j , ' n x (14) 

{Tre-P H °u(j3, 0)} 

To calculate G(t',j' | t, j), the evolution operator w(/3, 0) has to be expanded perturbatively, i.e., the expansion 
express of G(t',j'' | T,j) consists of terms, for instance, as 

—, Yl J iuh--- J in,j n drx--- dr„(f r [aUrOaj(T)aj i (Ti)a J1 (Ti)---aj n (r„)a >i (r„)]) , (15) 
n - ■ — ■ Jo Jo 

n ,31,— ,ln,3n 

with (O)o being average quantity with respect to the unperturbed part Hq of the Hamiltonian. Here, the Jij reads 

J, if i, j are nearest neighbors of each other, ,„ „. 

Jij = { (16) 
10, otherwise. 

Thus, the following quantities need to be calculated first 

G n 0) ( r i,*i;--- T ™,4 I n,ii;-"T n ,i n ) = (TV[at, (r{)a n (ri) • • • a t t «)a ln (t„)])o, (17) 

these are the n- particle Green's function with respect to Hq. Unfortunately, the Wick's theorem cannot be applied 
to calculate the above quantity, since the eigenstates of the non-perturbed part of the Hamiltonian ([3]) is local, and 
each annihilation operator in the expression should be paired with a creation operator on the same site, otherwise the 
result would be zero. Therefore, we should turn to the theory of linked-cluster expansion [H, [3, i.e. to expand the 

n-particle Green's function in terms of the cumulants Cm' '(r{, ■ • • , r' m \ t%, ■ ■ ■ r m ) = (T r [ a ^( T i) a ( T i) ■ 1 ■ a U T m) a ( T m)])o 
in which the particle operators are all on the same site. 
By defining the generation function as 

C [ ° ] [j,f] = ln(f T e ( ^ dr[r(r)a(T)+j(r)at(T)1) ) , (18) 
all the cumulants can then be calculated by 

2m 

(19) 



Cl 0) (ri , • ■ • , r' m | n, • • • r m ) - W) _ _ _ 5]{T , m)5r{Ti) . . . Sr{Tm) \j,f] 
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These cumulants can be used to decompose the above mentioned n-particle Green's function. As an example, the 
decomposition of one and two particle Green's functions is shown in the following 

GfV'y \T,i) = 5 i ,, i C[°\r' \t) (20) 
G 2 0) ( T i'*i; r 2^2 \ n,h;n,k) = S^.^s^^s^^c^i^,^ \n,T 2 ) 

+S^A^c[ 0) (r[\T 1 )C[°\^\r 2 ) 

+ < 5. ii , J2 ^, 1 C| 0) (rnr 2 )^ 0) (^|r 1 ) (21) 
In order to reduce the complication of the calculation, these cumulants can be represented diagrammatically as 

C[°\t>\t)= » » » (22) 

TIT 



Ct\<A I ti,t 2 ) = (23) 
meanwhile, the hopping parameter is represented diagrammatically as 



Jij = , » . ■ (24) 

i J 

With the help of the cumulants and the diagram rules mentioned above, each term in the expansion expression 
of G(r',j' | r,j) can then be represented diagrammatically easily. When considering the cancellation effect of the 
denominator in Eq. (fT4"]) . the one particle Green's function G(r',j' | r,j) only consists of connected diagrams. 

In perturbative calculation, the choice of the diagrams of the Green's function is very subtle. In the problem of 
investigating the phenomena closely related to the SF-MI phase transitions, choosing terms based on the order of J/U 
is not a good approximation, since in the regime of supernuid, this quantity is no longer a small quantity jl9j |. more 
importantly, the Green's function should be diverging around the phase transition point, however it is impossible to 
obtain such diverging behavior from a finite-order perturbation calculation, since it yields only a polynomial of J/U. 
Instead, we consider a resummed Green's function[13| which contains only the diagrams of single chains, in terms of 
Matsubara frequency, it is expressed as 

G 1 (u m ,i,j)= » m — »— + 1 (25) 

The reason for such a choice is following. Let us compare these two terms 



(26) 



although both of them are second order terms of hopping parameter, however, according to the topology of the 
underlying lattice structure, the former is smaller than the latter by factor of ^, d is the dimension of the system. In 
general, if a diagram contains n loops, it will be smaller than the corresponding single-chain diagram at least by factor 
of ^pr. Similar argument was stated in the case of fermions fl^ . In the limit of d — > oo, comparing to single-chain 
diagrams, all diagrams with loops may be neglected. In this sense, our choice of the resummed Green's function can 
be looked upon as a sort of mean-field treatment. However, from the above discussion, we see that, in principle, the 
Green's function method can easily be extended to a regime beyond mean-field. 

When translating this diagram expression to the expression in terms of cumulants and hopping parameters, through 
the Fourier transformation, we get 

Gi( Wm ,k) = (my (27) 

1=0 

where J(k) is Fourier transformation of the 

J(k) =J2 ^e lk rs e- lk rj (28) 
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and 



C[°\u; m ) = / C[ v >(r)e^dr 



(0), 



In the case of the triangle lattice, J(k) reads 

J(k) = 2 J 

Since 



cos^a) + cos(k x - + K y — - — J + cos^ — k x -) 



(29) 



(30) 



Cf (r) = (f T [at(r)a(0)]) o , (31) 
together with Z*- ' = Tr[e _ ^ H °] and Ho\n) = E n \n) (\n) stands for a state with n particles on the site), we then have 



In terms of Matsubara frequency, it reads 



n=0 



n + 1 



•En+1 — E n + i0J m E n — E n -1 + lOJri 



(32) 



(33) 



Since the temperature in experiment is extremely low, thermal fluctuations are negligible (l9j . In the limit of T — > 0, 
the system of unperturbed part falls into the ground state (suppose the occupation number in the ground state being 
n), and the momentum space Green's Function is then reduced to 



oo ^oo r 



1 f 

Gi(r'|0,k) = - Y,j_ 



1=0 

oc 



n + l 



-i i+i 



1=0 J ~ co v 7 



E n +i — E„ + icj m E n — E n _x + iuij, 
i 



• 1 — /j, + iuj n 



where /i 



where 



(n — n + iuj m )(n — 1 - M + «^ m )_ 
^, cD m = t' = r'{7. After a laborious yet straightforward calculation, we have 

r 4.0 ;=o v ^ y 



_ ' (Z + l)!(Z + fc)! 



(34) 



(35) 



(36) 



and n is the particle number on a site. Thus, together with Eq. (|30[) , the momentum distribution function rt(k) is 
reduced to 



oc , ,i 
l (k)=iV s mk)| 2 ^(- S l 
i=o ^ ' 



, , , „ . , a , V^i \ „ / , \/3a , o, . 
2cos(fc s a) + 2cos(fc x - + k y —^-) + 2cos(fc y — fc x -) 



(37) 



This is our analytical result. In the following, we are going to calculate the time-of- flight absorbtion picture and the 
visibility of the interference pattern. 
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IV. THE TIME-OF-FLIGHT PICTURES 

Recently, Becker et al. [H performed an experiment to create an optical triangular lattice and loaded the ultra-cold 
87 Rb atoms in it. The wavelength of the laser beams which they used to create the lattice is A = 830 nm, the 
corresponding lattice constant is a = |A. In fact, in the experiment, what they created was not a real 2D triangular 
lattice but stacked layers of triangular lattices, this stacking structure is formed by an extra standing wave laser beam. 
To eliminate the influence of the third dimension, they set the lattice depth in the third dimension to be Vq — 30Er, 
this setup makes the tunnelling in the third dimension negligible, and the system exhibits the property of 2D. In their 
experiment, they filled about 40 layers and the occupation for each layer amounts to about 4000, together with the 
chemical potential being /i = 133 nk, a good estimation of the occupation number per site is n = 2. 

In order to compare our analytical result with their experimental observation, the third dimension has to be taken 
into account. Due to the orthonormality of the Wannier function, the third dimension would not affect the in-plane 
hopping parameter, J can be calculated explicitly from the lattice potential ([7]) and the Wannier function (|5J| via 
Eq.flS), it reads 



J = Er 



-A- 

4 \2V 



V exv>- — (2V )i }. (38) 



However, the scattering behavior of the system would be 3D, i.e. in calculating U, the interaction strength should be 
taken as g = A.-Kh 2 a s jm (a s = 5.34 nm is the 3D s-wave scattering length of 87 Rb and m is the corresponding atomic 
mass) and the integral of the third dimensional Wannier function has to be performed accordingly. According to the 
experiment, by taking the lattice depth Vq = 30Er in the third dimension, a detailed calculation leads to 



U = 4£r^(7ry ) 1/2 (30) 1 / 4 . (39) 



With (|38|) and (|39|) as well as the Fourier transformation of the Wannier function (j8J) in hand, we plot the analytical 
expression Eq. ()37|) of the momentum distribution function n(k) for various Vq in Fig. [TJ Qualitatively, it is in a good 
agreement with what observed in the experiment [4J . 

2.62 2.7 2.82 2.92 4.0 



FIG. 1: The single chain Green's calculation of time-of-flight absorbtion picture for an ultra-cold Bose system in a triangular 
optical lattice for various Vo- 



V. THE VISIBILITY 



However, the time-of-flight absorbtion picture can only provide us qualitative impression of the SF-MI phase 
transition. In addition, due to various experimental conditions, for instance the sensitivity of the detectors in the 
experimental setup, it is hard to compare the theoretical result of the time-of-flight pictures directly to experimental 
ones quantitatively In order to winkle quantitative information out of the time-of-flight pictures, the so-called visibility 
[20l [21j has to be calculated. The visibility is defined as 

^max ^min / tr .\ 

v= , (40) 

From Eq. ([37l) . it is easy to find out that n max takes place at (0,47r/(\/3a)) and (2n/a, 27r/(V3a)) and other four 
equivalent points in the momentum space while n m i n is at (27r/ (y3a), 2it /a) , (v3a), 0), and other four equivalent 
points. All these points have the same distance from the original point, hence at these points the Wannier function 
w(k) takes the same value, thus the visibility is solely determined by the correlations at these points. The theoretical 
result as well as the experimental data taken from Ref.[4| are plotted in FigJ^i. 
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FIG. 2: (a) The visibility of the interference pattern against U/(6J). The solid blue line is the theoretical result from Eg. ((40 
and the red dots are from the experimental observation [J. (b) The dependence of U/(6J) on Vb- 



From FigJJk, it is easy to see that our analytical result shows that the visibility of the interference peaks of the 
system drops dramatically fast in the region of 10 < Vq < 12. In Fig|2p, We plot the dependance of U/(6J) on V , 
the corresponding range of Vq is from 2.65 to 2.82, as can be seen in the inset of FigJ^jD. The experimental estimation 
is that the corresponding region is from 3.5 to 4. There is roughly a 20% deviation between our analytical result and 
the experimental data, this is understandable: first, our calculation is a sort of mean-field result; second, during the 
calculation of J and U, the Gaussian approximation of the Wannier function is adopted. Nevertheless, our theoretical 
prediction is in a reasonable agreement with the experimental data. Moreover, from the figure of visibility, we see that 
the visibility is non-zero far deep inside the Mott state, this is due to the short-range coherence[21[ in Mott insulator. 

VI. DISCUSSION 

In conclusion, by treating the hopping parameter in Bose-Hubbard model as a perturbation, with the help of the 
re-summed Green's function method and cumulants expansion, the momentum distribution function of the ultra-cold 
Bose system in triangular optical lattice is calculated analytically. By utilizing it, the time-of-fiight absorption picture 
is plotted and the corresponding visibility is determined. The comparison between our analytical results and the 
experimental data from Ref.[4J exhibits a qualitative agreement. 

As we have seen in the discussion of the paper, this systematic approach can in principle be extended to regime 
beyond mean-field theory by adding loop diagrams, more precise expressions of Wannier function may be adopted to 
improve the accuracy of the analytical result. Moreover, the method presented here can not only be used to investigate 
homogeneous cold atomic system, but also can be used to systems with nearest neighbor repulsive interactions. Ultra- 
cold Bose gas on a triangular optical lattice accompanied by nearest neighbor repulsive interaction would exhibit 
geometrical frustration effect, hence our present work may shed some new light in this field. 
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